NONLINEAR RESONANCES IN THE SOLAR SYSTEM 



Renu Malhotra 
Lunar and Planetary Institute 
3600 Bay Area Boulevard, Houston, TX 
Tel: (713)486-2114; Fax: (713)486-2162 
E-mail: renu@lpi.jsc.nasa.gov 



Abstract 

Orbital resonances are ubiquitous in the Solar system. They play a decisive role in the 
long term dynamics, and in some cases the physical evolution, of the planets and of their 
natural satellites, as well as the evolution of small bodies (including dust) in the planetary 
system. The few-body gravitational problem of hierarchical planetary-type systems allows 
for a complex range of dynamical timescales, from the fast orbital periods to the very slow 
orbit precession rates. The interaction of fast and slow degrees of freedom produces a rich 
diversity of resonance phenomena. Weak dissipative effects — such as tides or radiation 
drag forces — also produce unexpectedly rich dynamical behaviors. This paper provides 
a mostly qualitative discussion of simple dynamical models for the commonly encountered 
orbital resonance phenomena in the Solar system. 
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1. Introduction 

Soon after Newton's formulation of the Universal Law of Gravitation it became evident that a dynamical 
model of the Solar System based upon a simple superposition of two-body motions was not adequate for 
the observations. The mutual interactions of the planets (and of their natural satellites) were necessary 
to fit the observations and infer various properties of the system. Theoretical efforts to determine the 
effects of planetary interactions led to the development of perturbation theory and averaging methods, 
and to many other developments in mathematics. In recent years, celestial mechanics has evolved from its 
traditional habitat in the study of orbits of particular Solar System objects to the study of the structure 
of extended regions of phase space (of planetary systems, of satellite systems, and of other small celestial 
bodies), as well as of evolution due to dissipative effects. The greatest difficulties in these studies — and 
also the most interesting dynamical behaviors — are associated with orbital resonances. 

Even in the simplest system consisting of only two planets, there are six degrees of freedom corre- 
sponding to the three spatial degrees of freedom for each planet. Therefore, what at first glance would 
appear to be only a two-frequency system (i.e., the frequencies of revolution of the two planets around 
the Sun) is actually one with six frequencies. Two of these are the obvious ones of revolution around the 
Sun; the other four (two for each planet) are the much slower frequencies of precession of the onentatton 
of the orbits. The existence of both fast and slow degrees of freedom produces a rich diversity of orbital 
resonance phenomena. 

As a general classification, there are two types of orbital resonances. The intuitively most obvious 
type, referred to simply as "orbital resonance" or "mean-motion resonance", occurs when the orbital 
periods of two planets - or satellites - are nearly commensurate. The second type, called "secular 
resonance" , involves commensurabilities between the slow frequencies of precession of the orientation of 
orbits. Classical perturbation theory runs into the notorious problem of "small divisors" in analyzing 
the mutual gravitational perturbations of two planets near a resonance. In the case of the first type of 
resonance, a significant but subtle difficulty is that there are actually several resonances in the vicinity 
of any orbital period commensurability. This "resonance splitting" arises due to the precession of the 
orientation of orbits. 

The multiplicity of resonances is a consequence of the coupling between perturbations of different 
timescales, and, in general, may produce chaotic behavior in the vicinity of an orbital resonance. However, 
in many cases, the behavior at a resonance is largely regular, albeit complicated by non-linear effects. 
Examples include the observed orbital resonances amongst the satellites of Jupiter and Saturn, where 
"single resonance theory" works [1], and many phenomena in planetary rings where resonant perturbations 
from satellites are implicated in a bewildering variety of features such as gaps, kinks and sharp edges, and 
the confinement of narrow rings [2, 3]. In other cases, interactions between neighboring resonances become 
important; these lead to a variety of different phenomena associated with secondary resonances and 
chaotic dynamics. A well-known example is the "Kirkwood gaps" in the asteroid belt, which are deficits 
in the number distribution of asteroids at several locations corresponding to mean motion resonances 
with Jupiter; chaos due to interacting resonances has been identified as a mechanism for producing 
at least one of these gaps [4, 5, 6] and for transporting meteorites from the main asteroid belt to the 
Earth [7, 8]. Chaotic resonances are of great significance in the formation of planetary systems [9, 10] - in 
planetesimal dynamics, the transport of planetesimals to a proto-planet, the interactions of proto-planets, 
and finally the clearing of planetesimal debris from interplanetary space [11, 12, 13]. An interesting and 
very complicated example is that of Pluto which exhibits several secular resonances in addition to a 3:2 
orbital period resonance with Neptune [14]. Secular resonances may play a determining role in the long 
term dynamical stability of a planetary system [15, 16]. 

Dissipative effects (such as gas drag and mass loss in the Solar Nebula, tides in the case of planetary 
satellites, collisions in planetary rings, and radiation forces in the case of dust particles) cause orbits to 
evolve across resonances, leading to separatrix-crossing phenomena and capture into resonance. Pluto's 
very peculiar orbit (it crosses the orbit of Neptune but is dynamically protected from close encounters by 
the resonance phase locking) may owe its origin to orbit evolution of the outer planets due to mass loss in 
the late stages of planet formation [17]. The preponderance of orbital resonance locks amongst the Jovian 
and Saturnian satellites is thought to be due to the differential orbit evolution induced by tidal dissipation 
[1, 18, 19, 20, 21, 22]. These satellites probably did not form in resonant orbits; however, once established, 
the orbital phase locking can be maintained for long times. Similar evolution probably occurred in 
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the Uranian satellite system also. However, owing to a small but crucial difference in parameters, the 
orbital resonances in this system were not long-lived [23, 24]. Uranus is less oblate than Jupiter and 
Saturn; in addition, the relative masses of the satellites and their orbital radii conspire to cause an 
insufficient splitting of resonances; as a result, the interactions between neighboring resonances give rise 
to an instability of resonant orbits. Tidal evolution and passage through resonance can have a profound 
effect on the geophysical evolution of satellites [20, 22, 23, 24, 25, 26, 27]. Orbital resonance phenomena 
are also found in the dynamics of circumplanetary and interplanetary dust particles that are perturbed 
by radiation and electromagnetic forces. For example, the dusty ring around Jupiter exhibits features 
that are attributed to "Lorentz resonances" (commensurabilities between the Keplerian orbit period and 
the spin period of the planetary magnetic field) that dominate the orbital evolution of charged dust grains 
[28]. As a final example, I can mention the recent prediction and detection of a "Solar ring" of asteroidal 
dust particles that spiral inward due to Poynting-Robertson light drag, and are captured into long-lived 
orbital resonances with the Earth [29, 30]. 

The diversity of orbital resonance dynamics precludes a complete treatment in these few short pages. 
Therefore, this paper is confined to giving a brief review of dynamical models for the simplest of the 
commonly encountered orbital resonance problems in the Solar System. No attempt has been made to 
give exhaustive citations to the literature, but the selected references should provide the interested reader 
with sufficient leads to the technical literature. 



2. Basic formalism 

A peculiarity of the unperturbed two-body bounded motion is its complete degeneracy (see, for example, 
[31]): the motion is periodic in all three spatial degrees of freedom, and the three frequencies are all 
equal (to the orbital frequency). This property allows a choice of several different coordinates in which 
the motion is separable. The variables that have been used traditionally to characterize the Keplerian 
motion are the classical orbital elements (see Figure 1)^: 

a semimajor axis 
e eccentricity 
i inclination 

SI longitude of ascending node 
w longitude of periapse 
A mean longitude 

These elements form a complete set equivalent to the six Cartesian components of the relative position 
and velocity vectors, r = Ri — Rq and v = Vi — Vo, in the two-body system. The coordinates and 
velocities can be transformed by means of a canonical transformation to action-angle variables (see, for 
example, section 143 in [33]); these are given as follows in terms of the orbital elements: 

Jl = fiVGMa, 01 = \ 

J'2 = Ji(l-\/l-e2), 92 = -zu (1) 

J3 = Ji(l - e2)i/2(l _ cosi), 6*3 = -fi 

where /j, = momi/A4 is the reduced mass, and A4 = mo -|- mi is the total mass of the two-body system 
consisting of the primary, mo and the satellite, mi.^ For a system with N > I satellites about the 
primary, one can define N sets of canonical variables in a coordinate system found by Jacobi (see, for 
example, p. 539-540 in [34]). These coordinates are mathematically equivalent to the set (J,0) for a 
hierarchy of N two-body systems (but of course, the ( J, 0) are no longer action-angle variables when the 

The definitions of and A may appear counter-intuitive, but arise naturally in celestial mechanics calculations. For a 
more complete description and motivation for this choice of variables, the reader may consult any celestial mechanics text, 
for example, chapter 6 in [32]. 

^In the context of the Solar System, the central body (either the Sun or a planet) is always much more massive than the 
secondaries (the planets or satellites, respectively). We use the term 'primary' for the massive central body and 'satellite' 
for the secondaries. 
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mutual interactions of the satellites are taken into account). The j-th two-body system consists of a total 
mass A4j = mo + mi + ... + nij and a reduced mass, /^j = mj(A4j — rrij)/ Aij, with relative position and 
velocity vectors, = Rj — Rein and Vj = Vj — Vc4i where (Rj, Vj) and (Rem Vc4i "'"'') are the 
position and velocity of rrij and of the center of mass of mo, mi, mj_i, respectively. 

The three-body problem, with a dominant central body and two smaller bodies in nearly circular, 
co-planar orbits, is a fundamental model for perturbation analysis in Solar System dynamical studies. 
The Hamiltonian for the planetary three-body problem, truncated to the second order in the masses, is 
given by 

H = Ho + Hu, (2) 
Gnionii Gmom2 
2^, 2^ 

Hi2 = -GmimA-. ^ r " (4) 

The unperturbed part, Hq, is the sum of a pair of two-body Hamiltonians, and H12 is the interaction 
potential. 

In the case of planetary satellite systems, an important additional perturbation is that due to the 
oblateness of the planet; the orbit-averaged oblateness perturbation potential is given by (cf. section 
11-3(C) in [35]) 



= - E ^^^^ {^) V - e|)-^/^(3 cos^ i, - 1), (5) 



where R and J2 are the equatorial radius and oblateness of the planet. 

The perturbation, H12, can be expanded in a Fourier cosine series in the angular variables, {\j , Wj , 
with coefficients that are power series in the usually small orbital eccentricities and inclinations. This 
expansion takes the following form [34, 36] 

H12 ~ - ^"^^"^^ Sk,q cos (f)k,q (6) 
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k ,q 



where 

ij>k,q = (k + q)\2 - k\i - (qizui + q2ZU2 + qs^i + q4^2), q = qi + q2 + qs + q4 (7) 
and, to leading order in e and i, 

5,,,^/,,,(a)el^^lel^=lil^^l4^^l (8) 

where a = 01/02 < 1 is the ratio of the semimajor axes of the two orbits, and fk,q{c() is a function of 
Laplace coefficients [34, 36]. The sum is taken over all integer values of k and qi. Note that each term 
of the series satisfies the so-called d'Alembert characteristic: the integer coefficient of or SI in the 
argument of the cosine corresponds to the power of e (resp. i) in the coefficient S. This property permits 
a choice of non-singular canonical variables for the analysis (see Section 3 and Eqn. 12 below). A further 
condition on the integer coefficients, q^-\- q^ = even, is also realized; this condition means that the lowest 
order terms involving the inclinations are of the second order. 

In practice, in most cases of interest only a few lowest order terms in the above series are used to 
analyze the dynamics. Of special note are the "secular" terms - those whose argument 4>k^q is independent 
of the mean longitudes, \j , and depends only upon the slow angles Wj and For orbital configurations 
that are far from low-order resonances, the secular terms alone (together with -ffob) provide an adequate 
description of the long-term behavior of the system: the semimajor axes remain unperturbed, while the 
shape and orientation of the orbits exhibit long-period variations, including a secular precession of the 
nodes and periapses. 

The secular motions of the nodes and periapses are important in separating the frequencies of several 
angular arguments 4>k^q having the same ratio of coefficients of the mean longitudes, Ai and A2. 

An orbital resonance exists if the relative orbital motion of two satellites is period. This happens 
when the orbital frequencies, rij = A j , are commensurate, i.e., the ratio n2/«i is close to a ratio of small 
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integers, p : p + q. The integer q is called the order of the resonance. In a resonant configuration, the 
longitude of the satellites at every q-th conjunction (i.e., when they pass each other) is nearly fixed or 
oscillates about a slowly varying center. Stable oscillation of successive conjunctions can be maintained 
only about certain well-defined longitudes. A qualitative feature of a stable resonance configuration is 
that conjunctions of the satellites occur away from the longitude at which the orbits are closest together. 
The possible centers of oscillation of conjunctions can be uniquely specified with respect to the periapses, 
Wj , or nodes, Slj . 

The periodic repetition of the geometry of the system enhances the effects of the mutual gravitational 
interaction of the satellites, so that over many successive conjunctions, their small mutual perturbations 
can accumulate into large (but usually bounded) changes of the orbital parameters. What follows below 
is a description of the conditions for exact resonance and the multiplicity of a p : p + q orbital resonance. 

Sufficiently close to a specific p : p+q resonance, the frequencies of some of the terms in Eqn. 6 (those 
with k = p) become very small. For g = 1 or g = 2, to second order in ei, 62, ii and i'j, there are six terms 
with very low frequencies associated with any orbital commensurability; each of these defines a unique 
resonance condition. (For larger values of q, there is an even greater multiplicity of the split resonances.) 
As illustration, consider the 2:1 and 3:1 commensurabilities. Table I lists the defining resonance condition, 
the corresponding resonance angle (fip^q (which oscillates when the resonance condition is satisfied, but 
rotates otherwise), and a nomenclature for these resonances. The latter is derived from the coefficient of 
the leading resonant term in the perturbing potential. The subscripts 1 and 2 in the table refer to the 
inner and outer orbits, respectively. 



Table I: Splitting of resonances 



resonance 


condition 


resonance 


angle 




'name ' 


of the resonance 


(4A2 


-2Ai 


- 2^2) = 


4A2 


- 2Ai 


-2^2 




4:2 


i'2 resonance 


(4A2 


- 2Ai 


- ^2 = 


4A2 


- 2Ai 


-^2 - 


- Qi 


4:2 


iii2 resonance 


(4A2 


-2Ai 


- 2fii) = 


4A2 


- 2Ai 


-2fii 




4:2 


i\ resonance 


(2A2 


-Ai - 


-wi)=0 


2A2 


- Ai - 


- Wl 




2:1 


ei resonance 


(4A2 


- 2Ai 


— W2 — 1^1) = 


4A2 


- 2Ai 


— W2 - 


— Wl 


4:2 


6162 resonance 


(2A2 


-Ai- 


- ^^2) = 


2A2 


-Ai- 


- ^2 




2:1 


62 resonance 


(3A2 


-Ai - 


- 2^2) = 


3A2 


- Ai - 


-2^2 




3:1 


i'2 resonance 


(3A2 


-Ai- 


- ^2 - = 


3A2 


-Ai- 


-^2 - 


Qi 


3:1 


iii2 resonance 


(3A2 


-Ai - 


- 2fii) = 


3A2 


- Ai - 


- 2fii 




3:1 


i\ resonance 


(3A2 


-Ai - 


- 2roi) = 


3A2 


- Ai - 


- 2wi 




3:1 


e{ resonance 


(3A2 


-Ai - 


- ro2 ~ "^l) = 


3A2 


- Ai - 


- W2 — 


Wl 


3:1 


6162 resonance 


(3A2 


-Ai - 


- 2^2) = 


3A2 


- Ai - 


- 2n72 




3:1 


62 resonance 



The form of the resonance angles suggests a linear transformation to "resonance variables" via the gen- 
erating function: 

T = Ii(^(p + q)X2 - pXi - q^i^ /q + l2{(p+ q)h - P><1 - q^2^/q 

+ h {(p + g)A2 - pAi - qwi^ /q + I4 (^(p + q)\2 - p\i - qw2^ jq 

+riAi + r2A2 (9) 

Using the definitions of the canonical variables in Eqns. 1, this transformation yields the following reso- 
nance variables: 

<i>\ = ({p+ q)\2 - pM - q^i^ /q, h = Hi\J GMiai{l - ef)(l - cosii) ~ ]^mi\/Gmoaii\ 
(f>2 = {{p + q)\2 - p\i - q^2j Iq, h = 1^2 \/gM 202(1 - e2)(l - cos 12) ~ ^m2\/Gmoa2il 



5 



(t>3 = (iP+q)^2-P^i-q^i)/q, h = /^i\/GA^iai(^l - \Jl - ef) ~ ]-^mi\/ GniQaiel 
<f>4 = ?)A2 -pAi - 5^2^/5, h = H2\/GM2a2{l - \J'^ - 4) - ]-^rn2\/ GniQa^el 

71 = Ai, Fi = //i^/GA^ifli + ^(7i +h+h + h) (10) 

72 = A2, r2 = ll2\/GM2a.2 - ^ (Ji + J2 + J3 + Ia) 

Averaging R over the non-resonant, "fast" variables, 71 and 72, one gets a Hamiltonian with 4 
degrees of freedom , (/)4); Fi and r2 are first integrals for the averaged Hamiltonian. Referring 
to Table I, it can be seen that the and (^2 degrees of freedom are coupled together by the 'mixed' 
ixi'i resonance; similarly for the (/)3 and (^^ degrees of freedom. Furthermore, to the second order in ej 
and i j , there are no resonant couplings across the eccentricity-type and the inclination-type resonances. 
Therefore, in a further simplification, one can consider a resonance Hamiltonian with only two degrees 
of freedom: for the inclination-type (or "vertical") resonances, or ((/)3, /a; (/)4, 74) for the 

eccentricity-type (or "horizontal") resonances. 

The resonances are wtll separated if the rates Wj , Slj are sufficiently different from each other that the 
resonance splitting is larger than the widths of neighboring resonances. In such cases, "single resonance 
theory" provides a good description of the resonance dynamics. In many other cases of interest, the 
degeneracy is only imperfectly removed, and the interactions amongst two or more neighboring resonances 
lead to secondary resonances and chaotic layers. These different types of behavior are discussed in the 
following sections. 



3. Single resonance theory 

In the simplest case of resonance, only one of the resonance angles is truly a slow variable. Then, one can 
once more average over the faster angles and set to constants the corresponding conjugate momenta; this 
reduces the number of degrees of freedom to one. The single resonance Hamiltonian has the general form, 
K(I, (f)) = Ko(I) + V(I, (f)). As the most common cases of interest in planetary and satellite dynamics 
involve small-to-modest orbital eccentricities and inclinations, it is sufficient to use a second-order Taylor 
series expansion of A'o(/), and retain only the leading term in V(I,(I)). The symmetry properties of the 
planetary three-body system yield the following single resonance Hamiltonian: 

A'(7, <j)) = vl- pi"^ + iii'liyl'^ cos q<j). (11) 

K{I, (f)) is a pendulum-like Hamiltonian. is a measure of the "distance" of the satellites from the exact 
resonance condition for unperturbed circular orbits; /3 measures the nonlinearity and /j, the strength of the 
resonant interaction between the satellites^. Figures 2 and 3 show the phase space in Poincare variables. 



21 cos (f) and y = V2Ism(f), (12) 

for the case q = I (first order resonance), and q = 2 (second order resonance). 

The topology of the phase space is as follows. All phase space orbits are periodic, but a critical 
trajectory, the separainx, whose period is unbounded, exists for i> greater than a critical value i>c given 

by 

3|//2;3/4|i/3 forg = l 

2\^\ for g = 2 ^ ' 

When it exists, the separatrix may have one or two branches; it divides the phase plane into two or three 
zones: an external zone with circulating orbits outside the outer branch of the separatrix; an internal 
zone inside the inner branch of the separatrix (when there are two branches); and a resonance zone 
between the two branches of the separatrix (when there are two branches) or inside the separatrix (when 



A slight notational degeneracy: the /i here should not be confused with the two-body reduced mass defined in the 
previous section. 
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there is only one branch). Except for a narrow range of parameters and initial conditions, most orbits in 
the resonance zone are Hbraimg oihits, i.e. the resonance angle, (f), executes finite amplitude oscillations, 
whereas most orbits in the external and internal zones are circulating orbits ((f) increases or decreases 
without bound). The half- width of the resonance zone is approximately, 

1/2 

(14) 

where 1^ ~ is the value at the center of the resonance zone; A7 defines the maximum deviation 

of the canonical momentum, I, from its exact-resonance value for which the resonance angle can still 
librate. 



A7 



2// 

J 



(2Ir) 



3/2 



Adiabatic evolution 

The principal effects of energy dissipation or external torques on satellite orbits are to evolve the orbital 
semimajor axes, and thus to make the parameter i> change with time. For sufficiently slow variation of i^, 
the evolution of the system (Eqn. 11) can be evaluated by means of the adiabatic theorem (e.g. section 
49 in [31]) and a few simple ideas regarding separatrix-crossing transitions. The reader may consult 
references [1, 37, 38, 39, 40] for detailed explanations; here we give only the highlights of the theory. 

The satellite system describes a trajectory that is coincident with a closed ("guiding") trajectory 
in the (x, y) phase plane of the frozen system with the instantaneous value of i^. As long as v does not 
change very much during the period of the guiding trajectory, the action is an adiabatic invariant; the 
action is proportional to the area enclosed by the guiding trajectory: 

J = ^ j I d^) = ^ j xdy (15) 

As v changes slowly, the phase space of the satellite system evolves. The guiding trajectory also evolves, 
slowly changing shape while preserving the enclosed area. The exception to the adiabatic invariance of J 
occurs when the guiding trajectory becomes nearly coincident with the separatrix. In this situation, the 
period of the guiding trajectory becomes very large, and hence the change in v in one period becomes 
increasingly large, violating the adiabatic assumption. As v continues to change, the guiding trajectory 
must cross the separatrix. During this transition, while the system is in the neighborhood of the sep- 
aratrix, an instantaneous action integral is not well defined. Nevertheless, in the subsequent evolution, 
the separatrix will evolve away from the guiding trajectory, and the new action will once again be an 
adiabatic invariant. 

In order to apply this theory to orbital resonances, it is useful to recall the following. 

- v (which measures the deviation from exact commensurability) is an increasing or decreasing func- 
tion of time according as the ratio of the orbital radii, a = ai/a2, is an increasing or decreasing 
function of time due to the external torques. 

- The areas enclosed in the resonance zone and in the internal zone (see Figs. 2 and 3) increase with 
i>. Also, the distance of the center of the resonance zone from the origin is an increasing function 
of i>. 



From the definitions (Eqns. 10), y2I is proportional to the orbital eccentricity (or inclination). 
Therefore, the average radius of a trajectory in the (x, y) phase plane is proportional to the average 
eccentricity (or inclination). 

For circulating trajectories (in the internal oi external zones), the action, J, is also proportional to 
the average eccentricity (or inclination). 

For trajectories in the resonance zone, J measures the libration amplitude; the distance of the center 
of the resonance zone from the origin is proportional to the average eccentricity (or inclination). 
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The adiabatic theory for passage through resonance can be applied directly to determine the evolution of 
the orbital parameters. As an example, consider the evolution across an ei-resonance. Let the initial value 
of the eccentricity be ei o and corresponding action Jq. Away from the separatrix, Jo is an adiabatic 
invariant as i> changes slowly. The values of the action on the internal and external branches of the 
separatrix are functions of i^, denoted by Jint(j^) and Jext(^), respectively. There are two cases to consider. 
In following the analysis below, it is useful to refer to the phase portraits in Figures 2 and 3. 

• v/vr decreasing from a large positive value (i.e. the satellites approach the resonance on 'diverging' 
orbits): Initially the resonance zone is located a large distance from the origin. Assuming that 
the guiding trajectory is initially close to the origin, it will lie in the internal zone whose total 
area decreases with time. At some value, v = i/rp ^ when the internal zone shrinks sufficiently that 
JintX^T) — Jo, the guiding trajectory makes a transition across the separatrix to the external zone. 
(A transition into the resonance zone is not possible as its area is decreasing with time.) The new 
value of the action, J^, is the value of the action on the outer branch of the separatrix at i> = i^rp, 
i.e. Jf ~ Jext(j^T)- Therefore, upon crossing the resonance, the eccentricity "jumps" to a larger 
value. 

• ly/iyr increasing from a large negative value (i.e. the satellites approach the resonance on 'converg- 
ing' orbits): In this case, initially there is no separatrix, and the guiding trajectory is a circulating 
trajectory. At i>c a separatrix appears with values of the action, Jint(^c) = 0, and JextX^c) > 0. If 
Jo < Jext(j^c), then the guiding trajectory evolves smoothly into the resonance zone. On the other 
hand, if Jo > Jext('^c), then a separatrix-crossing transition will occur at some value i>t greater than 
i>c when the guiding trajectory becomes coincident with the external branch of the separatrix. 

is given implicitly by the equation Jext(^'T) = Jo- In this case, transition may occur either into the 
resonance zone or into the internal zone because the areas enclosed by both zones are increasing 
with time. A probability of resonance capture can be calculated as a function of the initial action, 
Jo, or equivalently, as a function of ci o [38, 40]. 

In summary, one finds that capture into resonance (libration) can occur only when the satellites approach 
the resonance on converging orbits (increasing i>). Resonance capture is certain if the initial eccentricity 
is smaller than a critical value; for larger initial eccentricities, the probability of capture decreases. If 
capture does occur, then in the subsequent evolution the eccentricity increases secularly as the adiabatic 
invariance of the action "drags" the satellites along the resonance zone which moves away from the origin 
as i> increases. If capture does not occur, then the new trajectory will lie in the internal zone, with a 
final value of the action smaller than its initial value; i.e. the final eccentricity will be smaller than its 
initial value. 

The above analysis has been applied in modeling and inferring the tidal evolution of the satellites of 
Jupiter and Saturn [1], and also in a theory for the origin of Pluto's resonant orbit [17]. Other interesting 
applications are in the evolution of dust particles due to electromagnetic and radiation forces in various 
environments (the rings of the outer planets [28, 41], dust particles in circumsolar orbits[30]), and dusty 
circumstellar disks [42]), and the evolution of planetesimals and protoplanets in the gaseous environment 
of the early Solar system [43, 44, 45]. 

4. Interacting resonances 

We have seen that for any orbital period commensurability there are actually several possible resonances. 
In the best of circumstances, these resonances are well separated due to the secular motions of the apses 
and nodes, and one can analyze each resonance in isolation. The next level of complication arises from 
the interaction of two neighboring resonances. 

The simplest model for this interaction is obtained by treating the coupling between two resonances 
as a perturbation on the single resonance model (cf. Eqn. 11): 



The form of the perturbation term introduced here is suggested by the coupling terms in the interaction 
Hamiltonian, Eqn. 6 (see references [24, 46]). The frequency SI is approximately the frequency difference 




(16) 
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between the neighboring resonances. If SI is sufficiently large, the single resonance theory can be recovered 
by the averaging principle. However, as I increases, two properties of this system conspire to cause a 
breakdown of the single resonance theory: (i) the strength of the perturbation increases, and (ii) the 
strength of the primary resonance, and therefore its libration frequency, also increases, thus reducing the 
frequency gap between the neighboring resonances. This results in a chaotic broadening of the separatrix 
of the primary resonance, as well as the appearance of secondary resonances deep inside the resonance 
zone of the single resonance model. Secondary resonances occur at commensurabilities between the 
libration frequency and the rotation frequency of neighboring resonances. An example of this situation 
is shown in a surface of section in Figure 4. 

For initial conditions (I, <j)) inside the resonance libration zone, the Fourier decomposition of the 
perturbation term is given by 



where to is the unperturbed libration frequency of 4>; the amplitudes are exponentially small for large 
\k\. One can see from this that even when SI ^ w, close to the unperturbed separatrix, the secondary 
resonance condition, to = Q,/k, will be satisfied for sufficiently large k. In fact, for \k\ greater than some 
^min, all the secondary resonances will overlap and broaden the separatrix into a chaotic layer. The width 
of the chaotic separatrix is exponentially small with the ratio Q,/uj. 

When is not too much greater than to, low-order, isolated secondary resonances can appear near the 
center of the resonance zone (cf. Fig. 4). Single resonance theory can be constructed for these secondary 
resonances, with Hamiltonians of the same form as Eqn. 11 [46, 47]. 

Adiabatic evolution 

The appearance of secondary resonances near the center of the primary resonance zone leads to a new 
phenomenon in the tidal evolution of satellite resonances: an orbital resonance between two satellites can 
be disrupted by means of capture into a secondary resonance. This mechanism is described below. 

Resonance capture at a small value of (7) is described well by the single resonance theory as the 
chaotic layer near the separatrix is initially exceedingly small. Upon capture into the resonance, (7) 
increases while the libration amplitude becomes quite small. At several points in the subsequent evolution, 
as (7) continues to increase, the system evolves across several secondary resonances. These resonances 
are born at the center of the primary resonance and migrate out toward the separatrix as (7) increases. 
Capture into any one of secondary resonances drags the satellites towards the chaotic separatrix and 
eventually allows an escape from the primary resonance. A schematic diagram of this mechanism is 
shown in Figure 5. 

An example of such interacting resonances occurs in the Uranian satellite system. At present these 
satellites do not exhibit any orbital resonances, but temporary resonance capture in the past may help 
explain their inferred thermal history. Furthermore, it has been shown that a temporary 3:1 resonance 
between Miranda and Umbriel may account for the anomalously high inclination of Miranda [23, 24]. 
Figure 6 shows a numerical simulation of the tidal evolution through the Miranda-Umbriel 3:1 inclination- 
type resonances in which the satellites escape from the resonance by the secondary-resonance-capture 
mechanism. Miranda's temporary residence in a similar eccentricity-type resonance in the past may also 
help account for its geologically diverse surface, because a high orbital eccentricity would enhance the 
tidal heating within this icy satellite. 

5. Chaotic resonances 

The interaction of two or more resonances at the same orbital period commensurability produces chaotic 
layers at the separatrices of the individual resonances. When the frequency separation of the primary 
resonances is modest, the origin of the chaos is readily interpreted as the accumulation of secondary 
resonances near the separatrix (cf. Figs. 4 and 5). 



+ oo 




(17) 
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An instability of circular orbits near the 3:1 Jovian resonance in the asteroid belt arises due to the 
interaction of the ef and 6162 resonances, where ei is the eccentricity of an asteroid and 62 is Jupiter's 
eccentricity [4, 5, 6, 48]. Many previous investigations that assumed a circular orbit of Jupiter had failed 
to find any instabilities because Jupiter's eccentricity is essential to the origin of the chaotic zone at 
this resonance. The depletion of asteroids at some of the Jovian resonance locations is attributed to the 
chaotic excitation of high orbital eccentricity. Asteroids in highly eccentric orbits become Mars- or even 
Earth-crossers. Close encounters with these planets then clear out the 3:1 Kirkwood gap. It is natural to 
conjecture that this mechanism — chaotic dynamics in the planar, elliptic, restricted three-body system 
of the Sun-|-Jupiter-|-asteroid, together with close encounters with the terrestrial planets — may provide 
a universal explanation for all the Kirkwood gaps. The hopes for this have not been fully realized as 
yet because of significant differences in the phase space structure at other resonances (notably at the 2:1 
Jovian resonance) [48, 49]; a complete solution awaits future efforts. In order to solve this problem, it 
may be necessary to consider dynamical models more complex than that based upon the planar, elliptic, 
restricted three-body problem [6, 48, 49]. 

A large chaotic zone near the separatrix of a resonance may occur even in the case where the 
perturbing planet is on a circular orbit (i.e., in the planar, circular, restricted three-body problem) if the 
resonant orbit is sufficiently close to that of the planet. The origin of this chaos lies in the interactions 
between distinct mean motion resonances, for example, a 3:2 resonance and a 4:3 resonance. An example 
of this phenomenon is seen in a simple model for the Neptune-Pluto 3:2 resonance [14]. Pluto is modeled 
as a massless test particle (its mass is only two-thousandths that of Neptune), and Neptune's orbit is 
assumed circular (Neptune's eccentricity is currently only ~0.01, while Pluto's is 0.25). Pluto's inclination 
is also neglected in this model. It is found that the region of stable resonance librations has an amplitude 
of about 100 degrees and is surrounded by a chaotic zone (Figure 7). The observed libration amplitude 
of Pluto's known orbit is about 80 degrees, safely inside the stable zone. 

Evolution due to weak dissipation across chaotic resonances is a pressing and largely unsolved prob- 
lem. A handful of numerical studies on the tidal evolution of satellites across chaotic layers may be found 
in the literature [23, 24, 27], but only a few analytical results are available [50]. It is unclear how to 
quantify adtahattc evolution in this case, for a separation of the timescales of "chaotic diffusion" and of 
dissipation becomes difficult. For the same reason, numerical modeling of these types of problems is also 
fraught with difficulties. 

6. Concluding remarks 

The Solar System as a dynamical system is a source of a very diverse and large number of dynamics 
problems. Many of these have parallels in other areas of physics. This paper gives only a very brief 
glimpse of the simplest examples of these problems. The analytical approach taken here describes some 
of the basic mechanisms that operate in the dynamical evolution of planetary systems. Entirely omitted 
in this paper is a discussion of numerical modeling in celestial mechanics which is a very active area 
of research; indeed it exemplifies in many ways the state-of-the-art in computational precision available 
today. The interested reader may consult [51] and references therein. The very large parameter space 
and very large range of dynamical timescales in the planetary few body problem (with or without weak 
dissipation) continue to challenge analytical techniques as well as the fastest numerical modeling hardware 
and software. 
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FIGURE CAPTIONS 
Figure 1: 

This shows the two-body Keplerian orbit in space and the definitions of the orbital elements. SI is the 
longitude of ascending node, to is called the argument of periapse, and / is the true longitude. In the set 
of 'classical' orbital elements, the longitude of periapse is defined as w = uj + and the mean longitude 
A = J ndt + zu, where n is the orbital frequency. For a circular orbit, J ndt = f, but for an elliptical 
orbit, it is an "averaged" value of the true longitude. 

Figure 2: 

The phase space in Poincare variables (x, y) for an isolated first-order resonance for (a) v < v^, (b) v = i/^, 
and (c) v > v^. The orbits shown in (d) are the same as those in (c), but are displayed in the (/,</)) 
plane; the pendulum-like structure is apparent in these variables. 

Figure 3: 

Same as Figure 2, but for a second-order resonance. 
Figure 4: 

This is a surface of section for the perturbed resonance Hamiltonian (Eqn. 16) for a weak perturbation 
of frequency near three times the small amplitude libration frequency of the primary resonance. The 
apparition of secondary resonances and the chaotic broadening of the separatrix is evident. 

Figure 5: 

A schematic picture of the adiabatic passage through an orbital resonance. The horizontal axis is a = 
Oil 02 < 1, the ratio of the semimajor axes of the two orbits. Resonance capture is possible if the satellites 
approach the resonance 'from below', i.e. with a increasing with time. Capture in the orbital resonance 
is followed by a secular increase in the mean value of the canonical momentum, I; capture in a secondary 
resonance "drags" the orbit out towards the separatrix, allowing escape from the orbital resonance. 

Figure 6: 

A numerical simulation of the tidal evolution of Miranda through a 3:1 inclination-type resonance with 
Umbriel. Miranda's inclination is pumped up to a high value before capture in a 3/1 secondary resonance 
disrupts the orbital commensurability. (Reproduced from [24].) 

Figure 7: 

Surface of section in the circular, planar, restricted three-body model for the Neptune-Pluto 3:2 resonance. 
The variables are (f) = 3Ap — 2Ajv — t^p, and J = ^ap(l — ^1 — Cp). (The subscripts P and N refer 
to Pluto and Neptune, resp.; op is measured in units of Neptune's orbital radius which is assumed to be 
constant.) (Reproduced from [14].) 
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